Asymmetric orbital-lattice interactions in ultra-thin correlated oxide films 
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Using resonant X-ray spectroscopies combined with density functional calculations, we find an 
asymmetric bi-axial strain-induced d-orbital response in ultra-thin films of the correlated metal 
LaNiOs which are not accessible in the bulk. The sign of the misfit strain governs the stability of an 
octahedral "breathing" distortion, which, in turn, produces an emergent charge-ordered ground state 
with an altered ligand-hole density and bond covalency. Control of this new mechanism opens a 
pathway to rational orbital engineering, providing a platform for artificially designed Mott materials. 



PACS numbers: 73.20.-r, 73.50.-h, 68.55.-a, 71.20.Be 

Heteroepitaxial synthesis is a powerful avenue to mod- 
ify orbital-lattice interactions in correlated materials with 
strong electron-electron interactions derived from transi- 
tion metals with open d-shell configurations [TH2] . Epitax- 
ial strain allows access to latent electronic functionalities 
and phases that do not exist in bulk equilibrium phase 
diagrams [3H8]. However, efforts to rationally control 
properties that are exceedingly sensitive to small per- 
turbations [9] through the orbital-lattice interaction are 
impeded by the poor understanding of how heteroepi- 
taxy imposes constraints on the orbital response [TQifT2] . 
Despite the recent progress in strain-induced orbital engi- 
neering, a crucial fundamental question remains: when a 
single electron occupies a doubly degenerate d-orbital in 
a cubic crystal field as in perovskites with Cu 2+ , Mn 3+ or 
low spin Ni 3+ cations, how does the substrate imposed epi- 
taxial constraints dictate the correlated orbital responses 
of ultrathin films? 

The exceptional strain control of the frontier atomic or- 
bitals relies on the susceptibility of the orbital occupations 
and their energy level splittings to biaxial strain-induced 
lattice deformations, i.e. orbital-lattice interactions. The 
conventional orbital engineering approach in perovskite- 
structured oxides is often rationalized as follows: coherent 
heteroepitaxy imposes a tetragonal distortion on the film's 
primitive unit cell, which then modifies the chemical bond 
lengths of the functional octahedral building blocks. The 
bond distortions in turn alter the crystal field symme- 
try and remove the cubic twofold e g electron degeneracy. 
Since each e g orbital state is of the same symmetry, it is 
generally anticipated that both tensile and compressive 
strains should symmetrically alter the d x 2_ y 2 and d 3z 2_ r 2 
orbital states. They are either lowered or raised relative to 



the strain- free band center of mass, which for finite filling, 
leads to an orbital polarization. This symmetric strain- 
induced orbital polarization (SIOP) concept is routinely 
used to rationalize the orbital responses of many complex 
oxide systems [12] [13] and even theoretically suggested to 
be a possible route to replicate high-T c cuprate physics in 
nickelates [14]. However, the strategies for tuning orbital 
ground states in ultrathin films devised from the suppos- 
edly symmetric strain response of the orbital occupations 
alone are violated more often than they apply: Doped 
manganite thin films exhibit d 3z 2_ r 2 orbital polarization 
regardless of the bi-axial strain sign [15] [16]. Ultrathin 
cuprate bilayers also show variable critical temperatures 
[T7] correlated with the sign of the interface lattice misfit 
[TBI, an d cobaltite films are either ferro- or diamagnetic 
depending on the strain state [T9] . 

In this Letter, we provide insight into why many ultra- 
thin correlated oxides violate the SIOP model by investi- 
gating the orbital-lattice interactions in 10 unit cell thick 
films of the correlated and orbit ally degenerate (t^g^g) 
metal LaNiOs (LNO). This representative spin-^ system 
has no magnetic or structural transitions below 500 K, 
which otherwise might obscure the investigation [20H22] . 
We identify the microscopic mechanism responsible for 
the asymmetric strain-induced orbital response and show 
that it originates from latent instabilities of the bulk ma- 
terial, which are strongly enhanced due to the epitaxial 
constraints imposed by the heterointerface. We further 
demonstrate deterministic control of the ligand-hole den- 
sity and covalency through the asymmetric orbital-lattice 
response, enabling access to an emergent charge-ordered 
phase not attainable in the bulk. 

For this purpose, we synthesized high quality epitax- 
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ial LaNiOs films under different strain states by pulsed 
laser deposition on (001) -oriented single crystal LaAlOs 
(LAO, lattice mismatch of -1.08%), SrLaGa0 4 (SLGO, 
+0.31%), SrTi0 3 (STO, +1.83%), and GdSc0 3 (GSO, 
+3.6%) substrates. Synchrotron-based X-ray diffraction 
confirms the LNO films are fully and homogeneously 
strained [23]. Consistent with previous studies on rhom- 
bohedral perovskites [24 , 25 , we find our films have lower 
(monoclinic) symmetry than that of the bulk. 

Figure [TJa) reveals a large anisotropy in the galvano- 
electric response with the sign of the strain. On the 
tensile side, the inverse Hall coefficient R^ 1 is small and 
independent of lattice mismatch indicating nearly zero 
strain-modulation in the carrier density, while for com- 
pressive strain an order of magnitude increase in R^ 1 
occurs. In addition, Fig. |TJb) shows the thermoelectric 
power (TEP) over a broad temperature range from 300 K 
down to 1.5 K for LNO films under different strain states. 
The thermopower in films on the tensile side lack any 
noticeable anomalies and follow the Mott relation [26] , 
though the magnitude is shifted relative to the bulk to- 
ward positive values. For compressive strain, the magni- 
tude is comparable to the bulk at T>80 K, but a different 
temperature dependence emerges at low T: namely, a 
large negative peak, reminiscent of phonon drag, occurs 
at 25 K. The asymmetric strain-induced TEP response 
indicates distinct differences in the electronic structure 
under tensile and compressive strain. In contrast, hydro- 
static pressure produces only modest enhancements in 
the TEP magnitude [Fig. [TJb)] , indicating compressive 
bi- axial strain is inequivalent to hydrostatic pressure. 

To explore the microscopic origins of the asymmetric 
strain-induced electronic states, we probe the local oc- 
cupied 3d orbital configurations and the oxygen- derived 
states using both X-ray absorption spectroscopy (XAS) 
and X-ray linear dichroism (XLD) measurements at the 
Ni L- and O K-edges. Figures [2ja)-(c) show the normal- 
ized room temperature XAS spectra at the Ni Ls^ edge 
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FIG. 1. (Color online) Asymmetric strain-induced re- 
ponses: (a) inverse Hall coefficient, R^ 1 and (b) temperature- 
dependent TEP. The carrier-density modulation in perovskite 
oxides is intimately linked with the the local structural of the 
NiOe octahedra (inset). The TEP in (b) reveals that epitaxial 
strain and hydrostatic (bulk LNO under 1 bar and 13) kbar 
pressure are inequivalent. Bulk data from Ref. [32 . 



for polarization oriented parallel to a — b vs. c axes under 
different epitaxial strain states. Here, we focus on two 
representative films under (on) compressive (LAO) and 
tensile (STO) strain states (substrates). Each polariza- 
tion averaged spectra is representative of octahedrally 
coordinated Ni 3+ in good agreement with the valence of 
bulk LNO. There are distinct differences, however, in the 
polarization dependence between films under compressive 
and tensile strain: The absorption for photon polariza- 
tion perpendicular to the Ni02 sheets is shifted ~0.1 eV 
lower in energy than the in-plane polarization (a — b) 
for LNO under compressive strain on LAO [Fig. |2ja)]. 
The dichroism indicates a small conduction band splitting 
Ae g between the Ni-derived d x i_ y i and d^ z i_ r i orbitals 
(Fig. [2] right) due to the distortion of the spherical charge 
distribution around the Ni sites by non-uniform Ni-0 
bond lengths. The sign of the XLD implies that the in- 
plane Ni-0 bonds are moderately compressed relative to 
the out-of-plane apical Ni-0 bonds along the c-axis in 
agreement with both the model of symmetric SIOP and 
the measured axial ratio for the films on LAO. 

For STO with the same substrate-film orientation, but 
strain of the opposite sign, the SIOP model predicts an 
inversion of the orbital polarization and the population of 
the d x 2_ y 2 orbital; however, consistent with our measured 
R^ 1 and TEP, the polarization dependence of the films 
under tensile strain is remarkably different. Fig. [2^b) 
shows the orbital dichroism is absent, suggesting that 
uniform Ni-0 bond lengths and a highly symmetric cubic 
crystal field (Ae g ~ 0) remain seemingly intact despite the 
large (~ 2%) tensile strain. These findings are surprising 
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FIG. 2. (Color online) Normalized absorption spectra at the Ni 
L-edge (fluorescence mode) with varying photon polarization 
(see legend) for films under (a) compressive and (b) tensile 
strain, on LAO and STO substrates, respectively, and (c) bulk 
(ceramic) LNO. (right) Schematic orbital level diagram for 
bulk LNO and the inferred SlOP-effect on the e q doublet. 
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as the in-plane lattice parameters are expanded by 2% 
and due to coherent film registry with substrate, the 
symmetric SIOP should have produced a large orbital 
polarization. 

To elucidate how the strain-stabilized film structure pro- 
duces an asymmetric orbital polarization, we perform den- 
sity functional calculations within the local-spin density 
approximation (LSDA) plus Hubbard U method following 
the details in Ref. [25] We carry out structural optimiza- 
tions for homoepitaxially strained nickelate films using 
our experimentally determined lattice parameters. We 
find that samples under compressive strain have an atomic 
structure characterized by Ni-O-Ni bond angles, which 
are larger within the epitaxial xy-p\ane than those along 
the growth (z-) direction. The structure also consists of 
a single Ni 3+ site with different in-plane and out-of-plane 
Ni-0 bond lengths that lift the formerly bulk Dsd Ni-site 
symmetry. 

Films under tensile strain exhibit rotated octahedra but 
with an additional intriguing breathing distortion imposed 
on the Ni06 octahedra [Fig. [TJinset)]. The breathing 
octahedral distortion induces a Ni-0 bond length dispro- 
portionate and splits the formerly bulk equivalent Ni 
sites into Ni(l) and Ni(2). They are located at the cen- 
ter of nearly ideal NiC>6 octahedra with average in-plane 
and out-of-plane bonds lengths of 2.00 A and 1.93 A, re- 
spectively. Interestingly, an analogous tendency to bond 
disproportionation leading to charge disproportionation 
of the type Ni^ 3±5 ^ + is present in the low-temperature 
monoclinic phases of small rare-earth nickelates [27j [28] , 
but not for bulk LNO. Note, the Ni06 octahedral breath- 
ing distortion is unstable for compressive strain — a single 
Ni 3+ site is preferred. 

We next compare the electronic ground state of LNO 
on LAO and STO with the bulk electronic structure. For 
films under compressive strain [Fig.^a)], we find a strain- 
induced crystal field splitting of the Ni 3d orbitals with an 
electronic structure that remains metallic, as in the bulk 
[Fig. [3jb)] 5 albeit with a strongly increased bandwidth of 
~0.20 eV. We find that the low-spin Ni 3+ electronic con- 
figuration is accurately described as itinerant and strongly 
mixed with an oxygen ligand hole, d 8 L. In sharp contrast, 
films under tensile strain on STO show a semiconducting 
gap of 0.10 eV with the bottom of the valence band edge 
shifted by ~0.5 eV over the compressively strained film 
[Fig. |3^c,d)]. We observe a modulation in the Ni-0 bond 
covalency through the dpa band splitting, which produces 
an energy gap between the antibonding e g orbitals, that 
is stabilized by the emergent strain-induced octahedral 
breathing distortion not found in the bulk. 

Figures |3jc) and (d) show that the valence band of LNO 
under tensile strain is primarily composed of the e g states 
from the more ionic Ni(l) cation with nearly one electron 
in each d x 2_ y 2 and d% z 2 _ r 2 atomic orbital (S = 1). The 
suppressed intensity at the valence band edge indicates 
that this orbital occupation is due to charge transfer from 
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FIG. 3. (Color online) The majority spin Ni 3d density of 
states (DOS) for LNO on (a), LAO substrate; for (b), a strain- 
free bulk crystal; and on (c) and (d), STO substrate, (e) The 
electron density for LNO on STO in the MO2 epitaxial plane 
[values integrated between —0.70 and eV (denoted as valence 
band, VB) and and 1.7 eV (denoted as conduction band, 
CB)] reveals the 3D checkerboard spatial distribution of the 
Ni e g states derived from the small and large octahedral units. 
[Darker (blue) areas indicate a higher value of the charge 
density.] 

the covalent Ni(2) site to the larger and more ionic Ni(l) 
site. We calculate the charge transfer S between Ni sites 
to be ~0.20e and homogeneously ordered throughout the 
film [Fig. [5Je)]. Because 5 is small and the low-spin Ni 3+ 
valence is rare, it is equally likely that a bond-centered 
double ligand hole L state, similar to the Zhang-Rice 
states in hole doped cuprates [29] , could be stabilized over 
the site-centered configuration. Nonetheless, our XLD 
data and first-principles calculations show that for LNO 
films under tension, the preferred electronic configuration 
is one without a net d-orbital polarization. The semi-ionic 
state results from a charge density modulation that occurs 
in the presence of orbital level splittings induced by tensile 
strain; in contrast, compressive strain favors a uniform 
orbital polarization of the Ni sites. 

To evaluate the degree of strain-modulated covalency, 
we use resonant XAS measurements at the oxygen K-edge 
to isolate the change in the ligand hole density of the 
d 8 L component of the Ni 3+ ground state. The oxygen 
spectrum around 528.5 eV corresponds to transitions from 
the core- level O Is state into the unoccupied and strongly 
hybridized O 2p - Ni 3d states [Fig.|4^a)]; the observed pre- 
peak intensity provides a direct measure of the degree of 
Ni-0 bonding. Unlike the Ni L-edge which is well aligned 
to the bulk LNO value, the oxygen pre-peak reveals a 
remarkably systematic and nearly linear dependence of 
the peak position with strain [Fig.^b)]: it shifts at a rate 
of almost 100 meV per percent lattice misfit, indicating 
direct strain control of the ligand hole density. Fig. [4^c) 
shows a corresponding systematic change in the full-width 
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FIG. 4. (Color online) Bi-axial strain-induced covalency 
changes: (a) Low energy O K-edge pre-peak absorption around 
528.5 eV for compressive (red) and a range of tensile (shades 
of blue) strain states (a.u., arbitrary units). We omit peaks 
above 530 eV (transitions into La 5d4/ and Ni 4s4p states) 
for clarity, (b) Energy shift and (c) the normalized FWHM 
relative to the bulk LNO energy position with misfit strain. 

at half- maximum (FWHM), suggesting a strong narrowing 
of the electron bandwidth due to the deviation of the Ni- 
O-Ni bond angle away from the ideal 180° with increasing 
tensile strain. 

In charge transfer oxides like nickelates and cuprates the 
position of the O K-edge pre-peak directly couples to the 
charge transfer energy, which is also highly susceptible to 
the symmetry of the cation coordination [30]. The strain- 
stabilized atomic distortions should therefore modify the 
Madelung energy. The constant nickel valence across all 
strain states explored allows us to probe the limits of the 
cohesive energy tunability with X-ray photoemission spec- 
troscopy (XPS). Indeed, we find a controllable modulation 
of the O Is core level energy by ~ 0.29 eV between LNO 
on LAO and STO (0.16 eV per percent strain) due to the 
asymmetric orbital-lattice interactions in agreement with 
our XAS results and calculations [23] . These experimental 
findings demonstrate that strain control of the ligand hole 
density is due to a major Madelung energy response. We 
conjecture that by reversibly accessing these asymmetric 
orbital-lattice interactions control of the charge transfer 
energy and d-orbital polarization should enable access to 
latent metal-insulator phase transitions in other classes 
of correlated electron systems [3T] . 

In summary, using a combination of resonant X-ray 
spectroscopies and first-principles calculations, we report 
the discovery of a fundamental asymmetry in the het- 
eroepitaxial strain-induced orbital occupation response 
of ultra-thin perovskite films with orbitally degenerate 
transition metal ions. We suggest that knowledge of 
this asymmetric orbital-lattice interaction is fundamental 
to the rational design of quantum materials with exotic 
correlated phases and enhanced critical temperatures. 
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